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Thermodynamics of a three-level maser was studied in the pioneering work of Scovil and Schulz- 
DuBois [Phys. Rev. Lett. 2, 262 (1959)]. In this work we consider the same three-level model, 
but treat both the matter and light quantum mechanically. Specifically, we analyze an extended 
(three-level) dissipative Jaynes-Cummings model (ED-JCM) within the framework of a quantum 
heat engine, using novel formulas for heat flux and power in bipartite systems introduced in our 
previous work [E. Boukobza and D. J. Tannor, PRA (in press)]. Amplification of the selected cavity 
mode occurs even in this simple model, as seen by a positive steady state power. However, initial 
, field coherence is lost, as seen by the decaying off-diagonal field density matrix elements, and by 

the Husimi-Kano Q function. We show that after an initial transient time the field's entropy rises 
linearly during the operation of the engine, which we attribute to the dissipative nature of the 
£SJ , evolution and not to matter-field entanglement. We show that the second law of thermodynamics is 

satisfied in two formulations (Clausius, Carnot) and that the efficiency of the ED-JCM heat engine 
agrees with that defined intuitively by Scovil and Schulz-DuBois. Finally, we compare the steady 
state heat flux and power of the fully quantum model with the semiclassical counterpart of the 
^— I . ED-JCM, and derive the engine efficiency formula of Scovil and Schulz-DuBois analytically from 

\^») ■ fundamental thermodynamic fluxes. 
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I. INTRODUCTION 



Thermodynamics of quantum-optical systems has intrigued scientists ever since masers and lasers were realized 
experimentally. Scovil and Schulz-DuBois [l[ analyzed a three-level maser in the framework of a heat engine. Based 
on a Boltzmann distribution of atomic populations, they gave an intuitive definition of the engine's efficiency, and 
showed it to be less than an or equal to the Carnot efficiency. Using the concept of negative temperature [2|, and 
motivated by Ramsey's Q work on 'spin temperature', Scovil and Schulz-DuBois [4j extended their analysis of three 
level systems to cases where the reservoirs' temperature is negative, and introduced the concept of negative efficiencies. 
^ c*j Alicki studied a generic open quantum system coupled to heat reservoirs, and under the influence of varying external 
i] conditions (such as a time dependent field) [5(. Alicki partitioned the energy of a quantum system into heat and 
. work using the time dependencies of the density and Hamiltonian operators. Based on Alicki's definitions for heat 
and work, Kosloff analyzed two coupled oscillators interacting with hot and cold thermal reservoirs in the framework 
of a heat engine, and showed that the engine's efficiency complies with the second law of thermodynamics [||. In 
^ i' later work, Geva and Kosloff studied a three-level amplifier coupled to two heat reservoirs Q @]. In their model the 
external field influences the dissipative terms, and the second law of thermodynamics is generally satisfied. 

This paper is to some extent a continuation of the studies discussed in the previous paragraph. In contrast with 
previous work, in our approach the matter and the radiation field are treated as a bipartite system that is fully 
quantized, as opposed to a forced unipartite system. This treatment of the working medium (the material system) 
and the work source (the radiation field) on an equal footing requires some new thermodynamic developments, that we 
adapt from Q. The general methodology is applied to an extended dissipative Jaynes-Cummings model (ED-JCM), 
which consists of a three-level material system coupled to two thermal heat baths and a quantized cavity mode. We 
show that this system provides a simple model of light amplification, which can then be analyzed using formulations 
of the first and second law of thermodynamics for bipartite systems. The heat flux and power calculated with this 
model lead to an engine efficiency that is in quantitative agreement with the efficiency formula intuitively defined 
by Scovil and Schulz-DuBois. A semiclassical counterpart of the ED-JCM equations is then presented and solved 
completely at steady state, giving the efficiency formula of Scovil and Schultz-DuBois analytically from fundamental 
thermodynamic fluxes. 

This paper is arranged in the following manner. Section II is a brief introduction to the thermodynamics of bipartite 
systems. In Section III we define the ED-JCM master equation. In Section IV we present numerical results for the 
ED-JCM model, showing that it acts as a simple model for a quantum amplifier. In Section V we discuss the entropic 
behavior of the full system and its individual components, its behavior at steady state and the role of entanglement. 
In Section VI we give a thermodynamical analysis of the ED-JCM. We formulate the first law of thermodynamics in 
two different ways. We then show that the second law of thermodynamics is satisfied in two formulations (Clausius, 
Carnot), and that the efficiency of the ED-JCM heat engine agrees with that defined intuitively by Scovil and 
Schulz-DuBois. In Section VII we compare the steady state heat flux and power of the fully quantum model with a 
semiclassical version of the ED-JCM, and derive the engine efficiency formula of Scovil and Schulz-DuBois analytically 
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from fundamental thermodynamic fluxes. Section VIII concludes. 

II. THERMODYNAMICS OF BIPARTITE SYSTEMS 

A bipartite system is described by a density matrix of a C m ® C n Hilbert space. The partial density matrix of one 
part is obtained by tracing over the other: 

PA(B) = Tr B(A) {pab}- (1) 
The entropy of a quantum system is given by the von Neumann entropy [Iol | : 

S = -k B Ti{plnp}. (2) 
The evolution of a bipartite system is given by the following master equation: 

PAB = Ch[PAB] + Cd[PAB], (3) 

where Lh[pAB] = Pab] is the Hamiltonian part of the Lindblad super operator, and Cd[pAB) is the dissipative 

part of the Lindblad super operator. The bipartite time independent Hamiltonian is given by: 

H = H A + H B + Vab, (4) 

where H a = Ha ® Is and Hb = 1a ® H B are the Hamiltonians of subsystems A and B, and Vab is the coupling 
term between them. Here and throughout the article, we use bold letters to signify operators that have a tensor 
product structure. 

Heat flux and power of the individual parts of the system are defined by Q: 

Qa(B) = T*{Cd[PAB]H A ( B )} (5) 

Pa(b) = — -t^{pab[Ha(b)i Vab]}- 
The energy flux of the full system is due only to the dissipative part of the Lindblad super operator: 

E AB = Tv{C d [p A B]H}. (6) 

III. THE ED-JCM MASTER EQUATION 

Consider a three-level system interacting resonantly with one quantized cavity mode and two thermal photonic 
reservoirs as depicted in Fig. [TJ The system is governed by the following master equation in the interaction picture: 
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FIG. 1: Three level system interacting with two heat reservoirs (hot and cold) and a quantized cavity mode. 

Pmf = £-h[Pmf] + £dc[Pmf] + £dB:[Pmf}- (7) 
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The letters in the subscripts have the following significance: m=matter, /=field, <i=dissipative, /i=Hamiltonian, 
C=cold, H—hot. The Hamiltonian part of the Liouvillian is given by: 

i 

£h[Pmf] = -^[Vmf,Pmf], (8) 

where 

V mf = A(ct 2 i ® a) + a\ x ® a) (9) 

is a resonant JCM type interaction Hamiltonian, A being the matter-field coupling constant. CdciPmf] and CdH[Pmf] 
are the dissipative cold and hot Lindblad super operators, respectively: 

Cdc[Pmf] = r '02{(n02 + l)([<TO2 p m f , <T^ 2 ] + [<T 2 ,p m f tJ 2 ] ) +n 02 ( [crl 2 p m f , O- 02 ] + [Oo 2 , pmfCTOl] ) } 

£dH[Pmf] = r 01 {(n 01 + i)([croi Pmf, ^Ol] + [''"Olj Pm/ tr ol])+ n 0l([ <T OlPrn./ ' ""Ol] + [°01' Pmf cr Ol])}t (10) 

where F02 and Loi are the Weiskopf-Wigner decay constant associated with the cold and hot reservoirs, respectively, 
and no2 and noi are the number of thermal photons in the cold and hot reservoirs, respectively. Note that direct 
dissipation occurs only through matter-reservoir coupling (the cold photonic reservoir couples levels |0) and |2), the 
hot photonic reservoir couples levels |0) and |1)), and is typically used to represent atomic decay in quantum optics 
. The matter creation and annihilation operators are in tensor product form Uij — ov, ® 1/ , and their matrix 
form is given by: 

1 \ / 1 

'-21 = | <> I CT01 = I <7 2 = 

000/ 




The reservoirs' temperature is given by: 



C{H) fe s ln(l/n 2(oi) +1)' 



where ojc(h) is the central frequency of the cold (hot) reservoir. The ED- JCM master equation (equation [7]) can be 
obtained by summing the Hamiltonian contribution and the two dissipative contributions. Alternatively, it can be 
derived for a three- level system with a break in symmetry using the weak coupling (to the reservoirs), Markovian, and 
Weiskopf-Wigner approximations in a similar fashion to the simple JCM with master equation with atomic damping 
which is derived in Appendix I. 

The Hamiltonian (energy operator) of the full matter-field system is given by: 

H = H m + H f + V mfl (12) 

where H m = H m <8> 1/; H m — ha and Hf = l m <E> Hf, Hf — huifa^a are the matter and field Hamiltonians, 
respectively, and a is given by: 

w 
wi 
0J 2 

Under matter-field resonance (co m = uj\ — u) 2 — tOf) the Hamiltonian in the interaction picture is unchanged and is 
not time dependent (H 1 = exp% Hot H exp"i Hot = H, H = H m + Hf) since [H Q , V m f] = (when there is no 
resonance one can still transform to an interaction picture in which the Hamiltonian is unchanged [l2j]). However, as 
indicated previously (eq. [5]), in the interaction picture, the Hamiltonian part of the evolution of the density matrix is 
only via the interaction term V m f = H 1 — Hq. 

Before we move on to discuss the ED-JCM as a quantum amplifier, we wish to discuss the main differences between 
the ED-JCM and the quantum theory of the laser due to Scully and Lamb (SL) [l3| [3| ■ Firstly, in the SL model 
the material system (the atom) has either four levels [l4| or five levels [TTj] , whereas in the ED-JCM the matter 
has three levels. Secondly, in the SL model the transitions between the two upper lasing levels and the two lower 
levels is achieved through a phcnomcnological decay, whereas in the ED-JCM population may also be pumped from 
the ground state to the two upper lasing levels through the full dissipative Lindblad super operator. Thirdly, in 
the SL model the atom is assumed to be injected into the cavity in the upper lasing level and interact with the 
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cavity for a time r, whereas in the ED-JCM the matter is in continuous contact with the quantized cavity mode, 
and amplification is achieved for a wide range of initial states. Finally, in the SL model the field is allowed to decay 
using the Weiskopf-Wigner formalism, whereas in the ED-JCM discussed in this paper the field does not decay. In 
principle, cavity losses can be introduced to the ED-JCM. However, we do not consider field damping in this paper, 
which allows us to compare the thermodynamical fluxes in the quantum ED-JCM with their analog in a semicjassical 
ED-JCM (section VII) and a similar model by Geva and Kosloff in which field damping is not included The 
differences between our model and that of the SL model will be seen below to play a crucial role in our ability to give 
a thermodynamic foundation of amplification. 

IV. THE ED-JCM AS A QUANTUM AMPLIFIER 

The ED-JCM master equation, eq. was solved using the standard Runge-Kutta method (fourth-order for 
various choices of parameters. The accuracy of the solution was checked by decreasing the step size. Furthermore, in 
order to test whether the numerical solution captures all time scales (especially the rapid oscillations), the algorithm 
was tested on the simple JCM [lj| which can be solved analytically [TtJ [HI- In all plots presented here To2 = Toi = 
r = 0.001, A = 1, no2 = 0.1, uqi = 10, and quantities are given in atomic units. The condition A ^> T corresponds 
physically to a situation where the coupling between the matter and the selected quantized cavity mode is much 
stronger than the matter-reservoir coupling. 

The energy flux of the full matter-field system and the individual subsystems is given by: 

E mf ee Tr{p mf H} = Tr{C d [ Pmf }H} 
E m EE Tr{p m H m } = -jTv{p m f[H m , V m f}} + Tr{C d [pm.f}H m } = P m + Q m 

E f ee Tr{p f H f } = ~Tr{ Pmf [H f , V mf }} = P f , (13) 

where E m f, E m , and Ef arc the energy fluxes of the full matter-field system, the matter, and the field, respectively. 
p m and pf are obtained from p m f by a partial trace over the field or a partial trace over the matter, respectively. H m 
and Hf are the Hamiltonians of the matter and field subsystems, respectively, without the tensor product with the 
identity. Note that the energy fluxes of the individual subsystems in eq. [13] are defined via p m and pf together with 
the subsystem Hamiltonians H m ,Hf. In the next two subsections we discuss the transient and steady state energetic 
behavior of the ED-JCM. Since there is direct dissipation only through matter-reservoir coupling, there is no heat 
flux associated with the field (this is physically expected, and was shown analytically elsewhere [9(). 

A. Transient behavior 

The energy of the full matter-field system and of the individual subsystems is plotted in Fig. [2] for an initial state 
where the matter is in state |1) and the selected cavity mode has no photons (p m f(0) = (|l)(l|) m (g> (|0)(0|)y). At 
short times, t < TJ^, the matter and field energies oscillate at a frequency of SI = ir\ [l7|. Here T e ff — r " 01 "^" 02 is 
the effective decay constant. 

Moreover, at short times the well known collapse and revival phenomena [li| [2(| is observed for a sufficiently 
excited initial coherent state as depicted in Fig. [3] 

In order to monitor the field's coherence we calculate the quantum optical Husimi-Kano Q function which is defined 
by HI: 

Q(a r ,ai) = -(a\p f \a), (14) 

where \a) is a (generally complex) coherent state. In Fig. 3] we plot the Q function at four different times (t = 0, 
t = 0.026r~ / 1 / , t = 0.4T-} f , t = 253r ( T / 1 / ) for the initial state p m f(0) = (|l)(l|) m ® (l«)H)/, M 2 = 5- The Q 
function at t = is that of a coherent state with the phase centered around the real axis (Fig. [4Jl). At transient 
times, the Q function spreads in phase space, but it is not homogeneous as seen in Fig. UJd and Fig. 

B. Steady state behavior 

At much longer times, t > TJ^, the matter energy decreases to a steady state value, while the field energy increases 
with a steady state power of Pf s = 4.5975e -5 (the numerical value of a linear fit to the last 10000 points, R 2 = 1.000) 
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FIG. 2: Energy of the full matter-field system (solid line) and the individual subsystem (field dotted line, matter dash-dot line) 
for an initial state where the matter is in state |1) and the selected cavity mode has no photons (p m /(0) = (|l}(l|) m ®(|0}(0|)/). 
Note that at long times there is a steady state increase in the field's energy. 
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FIG. 3: Matter energy at short times (t < (2r) 1 ). The revival time (top) is t r = 2,rm V <") . ^ _ | a |2 _ ^qq^ wnere m j s 
positive integer. The Rabi oscillation time (bottom) is iR a bi = 



as seen in Fig. [5] Another indication for an increase in the field's energy is seen in the steady state increase in the 
full system energy. Thus, the field and the matter-field system as a whole never reach a steady state for the type 
of evolution discussed in this paper. The relation between the full system energy and the energy of the individual 
component subsystems will be discussed in the next sections. A steady state increase in the field's energy is clearly 
an amplification of the selected cavity mode. This behavior contrasts with the simple JCM in which the atom and 
field oscillate forever (the atom oscillates between the excited and ground states while the field oscillates between the 
|0) and |1) Fock states). Amplification of the selected cavity mode will occur with any other coherent state, including 
the |0) Fock state. The fact that the field's energy increases monotonically is not unreasonable, since the harmonic 
oscillator is infinite, and since we do not consider direct dissipation of the cavity mode (which could be modeled by 
transmissive mirrors if desired). 

The collapse and revival phenomenon at longer times is completely damped due to the dissipative contribution to 
the Liouvillian as seen in Fig. [3J At these long times all phase (internal coherence) information is lost: the Q function 
is radially symmetric and is dispersed on a bigger area (bottom of Fig. |4j;). From this time onwards the shape of 
the Q function remains unchanged, and it expands fully symmetrically. The decay of the initial field coherence is 
also reflected in the decay of the off-diagonal field density matrix elements. At t — lOrj^-, the off-diagonal matrix 
elements are 10 -13 times smaller than their initial value, and are practically zero. All the remaining density matrix 




FIG. 4: Husimi-Kano Q function of the selected cavity mode, p/(0) = a)(aj; \a\ 2 — 5. a): at t — the Q function is a narrow 
2D gaussian. b) and c): at t = 0.026F~^, t = 0AT~^ the Q function is spread in phase space inhomogeneously. d): at 
t — 253r~j,j> the Q function has expanded (due to amplification of the selected mode) into a radially symmetric annulus (all 
the initial phase information is lost). 



elements are diagonal with a Poissonian-like photon distribution whose average number of photons increases with 
time. 

The full density matrix can be divided into a 3 x 3 block matrix, each block associated with one element of the 
matter density matrix. At long times, t > 2T~^p the matter-field inter-coherence is maintained by the non- vanishing 
matrix elements: pi ; n,2;n+i> />2;n+i,i;n of the full density matrix. These elements correspond to matter-field coupling, 
maintained via the structure of the JCM Hamiltonian. Other density matrix elements at these long times are 5 — 9 
orders of magnitude smaller than the dominant non vanishing elements discussed above. 



V. ENTROPY IN THE ED-JCM 



We consider now the entropy in the ED-JCM. In the next two subsections we discuss the transient and steady state 
entropic behavior. In Subsection C we discuss the relation between the entropies of the individual subsystems and 
entanglement, both at transient and steady state times. 



A. Transient behavior 



In Fig. [5] we plot the entropy of the full matter-field system and the individual subsystems for the initial state 
Pmf(0) = (|l)(l|) m ® (\a)(a\)f, \a\ 2 — 25. The entropy plots at the top of Fig. [5] show that there is a rapid rise 

(t k — |— ) in the entropies of the matter (dash-dot line), the field (dotted line), and the full matter-field system 
(solid line). This overall rise in entropy is discussed in Subsection C) and is attributed to the dissipative nature of 
the problem. 

At times t < , the entropy of the individual subsystems (matter or field) is oscillatory, as seen by the matter 
entropy plot in the middle of Fig. [5] This behavior is typical of the simple JCM [TtJ At this stage S m +Sf > S m f, 
in Subsection C) we attribute the excess entropy to entanglement. 
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FIG. 5: Entropy of the ED-JCM for the initial state p m /(0) = (|l)(l|) m ® (|a)(a|; |a| 2 = 25)/. a): entropy of the full matter- 
field system (solid line) and the individual subsystem (field dotted line, matter dash-dot line) for the full evolution, b): matter 

entropy for times t < . c): entropy of the full matter-field system for times t > 2T~^. 



B. Steady state 



Fig. suggests that at t > rrjy the matter-field system has reached a steady state. Indeed, at times t > T~f, the 
energy of the matter remains constant (dash-dot line in Fig. [2]). However, as was indicated in the previous section, 
the field energy plot (dotted line in Fig. [2|) and the matter-field energy plot (solid line in Fig. [2]) both show a constant 
rise for t > ^ff Furthermore, a closer inspection of the matter-field entropy (Fig. Hh) reveals a constant slight rise 
in entropy at times t > \ (a similar rise in the field entropy is also observed). Moreover, the field density matrix 

eigenvalues change in the second and third significant figures over a F~}f time scale. These findings give further proof 
of the fact that the field and matter-field system as a whole never reach a steady state. 



C. Entanglement 



We will now analyze the nature of the entropies associated with the subsystems. The entropy of the individual parts 
of a bipartite system is closely tied to the issue of entanglement [22| [H| • An important measure for entanglement is 
the conditional entropy, defined for the matter-field system by: 

S(m\f) = S m f - Sf 

S(f\m) = Smf — S m , (15) 

where S(m\f) is the conditional entropy of the matter, and S(f\m) is the conditional entropy of the field. In contrast 
with the conditional entropy in classical bipartite systems, the conditional entropy in quantum bipartite systems can 
assume negative values. In this case, the correlation between the two parts of the system is of a purely quantum 
nature, and the system is therefore entangled. A fine example for entanglement in the context of our work is the simple 
JCM. Consider an initial state given by: p m f(0) = (|e)(e|) ® (|0)(0|)/, where the atom (indicated by subscript a) 
is in the excited state and the cavity mode is empty. Under the JCM Hamiltonian (pure Hamiltonian dynamics), the 
full atomic-field entropy is constant S a f{t) — 0. However, during most of the evolution time the conditional entropies 

of both the atom and field (which are equal) are negative. In our case, at short times, < t < we find that the 
matter's conditional entropy is negative. Thus, the excess entropy S m + Sf > S m f at these times is attributed to 
entanglement. 
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A more powerful test for entanglement, introduced originally by Peres [23J], is the negativity of the partially trans- 
posed density matrix. The partially transposed density matrix is defined by: 

Pia,j/3 = PiP,i<*- (16) 

A sufficient condition for entanglement is the negativity of p T2 . However, since this test applies only to finite 
dimensional density matrices, one should take care not to mistake truly negative eigenvalues with negative eigenvalues 
that are an artifact of truncation of an infinite Hilbert space [l8| . Indeed, at times smaller than the typical decay time 
(r = r^yy), we find that the matter-field partially transposed density matrix p^f is negative (negative eigenvalues 

3r _1 

with a substantial absolute value are found up to a time t — f"-, and hence the matter-field system is entangled. 

At t > F~}f the matter and field conditional entropies are positive. Moreover, the conditional entropies are almost 
equal to the partial entropies: s (™\f> , g ^ m ) > 99%, and p^j, is positive (as was indicated before). All these findings 
lead us to conclude that in all likelihood at long times the matter-field system is only weakly classically correlated. 

We summarize this section by stating that at short times (t < T~^), when the partial entropies are oscillatory 
(see bottom of Fig. [5J, the matter- field system is entangled, as verified by the negative conditional entropies and the 
negative partially transposed full density matrix. However, as dissipation sets in, the matter-field system becomes less 
and less entangled. At t > r^".., when the partial entropies are not oscillating any more, the matter- field system in 
all likelihood is not entangled (as verified by the positive conditional entropies and the positive partially transposed 
full density matrix), and the overall rise in entropy is attributed to the dissipative Lindblad super operator. 

VI. THERMODYNAMIC ANALYSIS OF THE STEADY STATE SOLUTION 

A. The first law 

The first law of thermodynamics is essentially given in equation 1131 However, some fine details need more clarifi- 
cation. The first law of thermodynamics for the full matter-field system in differential form is given by: 

E mf = Tr{p mf H} = Tr{C d [p mf ]H} = Q m + Q f + Q v = Q m + Q v , (17) 

where Qf = Tr{Cd[p m f}Hf} = as was shown elsewhere [9J, and Qv = Tr{£<j[p m ^] V m f}. Q m = 
TT{Cd[p m f]H m } = Tr{C d c[Pmf]H m } + Tr{£ dH [p m f]H m } = Q mC + Q m H is the heat flux associated with the 
matter and it is composed of heat fluxes from/to the cold and hot heat reservoirs. Note that to an observer looking 
on the matter-field system as a whole, the full system is only dissipating heat. 

Another way to formulate the first law of thermodynamics is based on the energy flux of individual subsystems. 
The first law of thermodynamics for the matter and field separately (in differential form) is given by: 

E m = Tr{p m H m } = ~^Tr{p mf [H m , V mf ]} + Tr{C d [p mf ]H m } = P m + Q rn (18) 
E f ee Tr{p f H f } = ~Tr{p mf [H f , V mf }} = P f , (19) 

where P m = — j{YY{p m f\H m , V™,/]}, and Pf = —j-Tr{p m f[Hf, Vmf]} are the power terms. Since we are considering 
the case of perfect matter-field resonance, P m = —Pf ([H mi V m f] = — [Hf, V m f]), hence: 

E m + E f = E mf -Q v . (20) 

It may be shown that Qv vanishes if the off-diagonal matrix elements of p m f are purely imaginary. Note that to 
an observer looking on the matter alone work flux (power) and heat fluxes are identified according to eq. 1181 in 
agreement with the traditional thermodynamic partitioning of energy into work and heat. The field, which is the 
work source, either receives or emits energy to the working medium (the matter) in the form of power. In this paper 
we are interested in optical amplification. Under such conditions, at steady state the energy flux balance is such that 
< 0, QmH > 0, QmC < and the three-level system operates thermodynamically as a heat engine. 
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B. Second law. Clausius formulation 

The second law of thermodynamics is obtained via the entropy production function of the full bipartite matter-held 
system, which is defined by (24[, 

+ (21) 

where dS ^ f is the entropy production associated with the bipartite matter-field density matrix (via differentiation of 
the von Neumann entropy), and J is the entropy production associated with the reservoirs (via the heat flux from/to 
the reservoirs) given by: 

J = -PcQc - PhQh , (22) 

where (3 C (H) = fe^q/f))" 1 . and Qc(H) = ^{^dC(H)[Pmf](H m + V mf )} = Q mC (H) + Qvc(H)- Spohn showed that 
for a completely positive map (such as the Lindblad super operator) [241 ]: 

a Q > 0. (23) 

Equation 1231 represents the differential form of the second law of thermodynamics in Clausius's formulation, since the 
sum of the entropy changes of the system and reservoirs is guaranteed to be positive. 

C. Second law. Carnot's formulation 

We now define a new entropy production function: 

a m = -7^- + J m , (24) 

where is the entropy production associated with the matter density matrix (via differentiation of the matter von 
Neumann entropy), and J m is the entropy production associated with the reservoirs, 

Jm = —PcQmC — PhQiuH, (25) 

taking into account the contribution only from the matter heat flux 

Qm = QmC + QmH- (26) 

The physical idea behind a rn is that it is built only from matter thermodynamic fluxes: the intrinsic entropy flux 
^f- and the entropy flux J rn arising just from matter heat fluxes. For many initial matter states Qvc(H) = at all 
times, and hence J m = J. Moreover, when the matter reaches a steady state we always find numerically (irrespective 
of the initial matter state) that QvctH) — 0. This is the case also in the semiclassical ED-JCM discussed in section 
VII, where it can be shown analytically that QyciH) = ®- Therefore, at steady state, a m is physically similar to 
the entropy production function in the semiclassical case, a sc , where the field is not quantized. In contrast with 
a sc and cr^, a m is not guaranteed to be positive at all times (especially at times t < (2r) _1 , due to the highly 
oscillatory nature of the partial entropy at short times). However, when the matter reaches a steady state (^ff- = 0), 
the increase in the field's entropy is marginal (as was indicated before), and the main source of entropy production is 
the heat flux from/ to the heat reservoirs ( J > ^ f ) . Thus, when the matter reaches a steady state 

o m = J m = J > 0, (27) 

and since the matter operates in a heat engine mode (Q m H > and Q m c < 0), we obtain Carnot's efficiency formula: 

Pa QmC + QmH . Tr - Tc - . 

V = —X = a (28) 

QmH HmH 1 H 

where we have used eq. [T8]with eq. [26l and eq. [27] with eq. [25] For example, the efficiency of the heat engine for 
the choice of parameters discussed in the previous plots and for various initial field strengths (ranging from an empty 
cavity up to 100 photons) is 75%, which is less than the Carnot efficiency which is 99%. 
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Scovil and Schulz-DuBois gave an intuitive, but non-thermodynamics definition of the efficiency of the three-level 
system operating as a maser |l|: 

Vm = —, (29) 

Up 

where u s is the signal (maser) frequency, and uj p is the pump frequency (central frequency of the hot reservoir, 
uj x — u ). By substituting our initial choice of parameters (u s = u = 0.075, and = 0.1) we see that 

our numerical result agrees precisely with the efficiency estimated by Scovil and Schulz-DuBois. We find numerically 
that the efficiency calculated from eq. [25] agrees precisely with the efficiency calculated using eq. [59] Although we do 
not have analytic proof of this equivalence for the case of quantized light, we show below that this equivalence can be 
derived analytically when the light is treated classically. 



D. The ED-JCM: A work source with an entropy content 



A work source is the physical entity on which work is done, or which performs work on a system (working medium). 
Conventional wisdom in classical thermodynamics states that a work source's entropy is constant during the operation 
of a heat engine [25[ . Whether the classical engine operates cyclically, as in the usual Carnot cycle, or synchronously, 
the working medium returns to its initial state. The working assumption in thermodynamics is that entropy may be 
produced at the boundary of the working medium and the heat reservoirs, but not at the boundary with the work 
source. 

The work source in the quantum amplifier discussed in this paper is the selected cavity mode which is amplified. 
At steady state, the density matrix of the matter becomes constant and thus its entropy is unchanged from this time 
onwards. Since energy is flowing from the hot reservoir to the cold reservoir and work is produced in the form of 
amplification of the cavity mode, this corresponds to the engine operating in synchronous mode with the cavity mode 
as the work source. Inspection of Fig. [5] shows that the entropy of the light is not constant. Even after the matter 
reaches a steady state, the entropy of the light continues to grow linearly in time. 



VII. THE SEMICLASSICAL ED-JCM 



A. Equations of motion 



The semiclassical ED-JCM master equation is similar to the quantum ED-JCM master equation given by equation 
[7J However, since the selected quantized cavity mode is replaced by a time dependent field, major differences arise. 
The field is considered as an external degree of freedom, and hence it has no entropy content. We propagate a 3 x 3 
density matrix representing the matter only, and all operators are represented by 3 x 3 matrices (as opposed to 
(3<8>n) x (3(gm) in the fully quantized case). Finally, the Hamiltonian part of the Liouvillian assumes a different form, 
where the creation and annihilation field operators are replaced by clockwise and anti-clockwise oscillating exponents. 
Despite the last difference, we note that in perfect matter-field resonance the Hamiltonian part of the Liouvillian in 
the interaction picture is time independent. The semiclassical Hamiltonian is given by: 

H = H m + V, (30) 

where H m is the matter Hamiltonian as given in equation 1121 (without the tensor product with lj), and 

V = A sc (a 2 ie iwt + 4 ie - iu)t ) (31) 

is the interaction Hamiltonian with a classical single coherent mode in the RWA. X sc is the semiclassical matter-field 
coupling constant (which can be obtained via the semiclassical coupling matrix element (26j ) given by (atomic units): 

A sc = D • (32) 

where D is the dipole operator, e is the field polarization, and Eg is the field amplitude which can be estimated by 
calculating the average value of the quantum field operator for a coherent state [26j : 

£o=(^)" /2 M, (33) 
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where u> is the mode frequency (not necessarily in resonance with the atomic transition) , V is the cavity volume, and 
| a | is the field strength. The quantum matter-field coupling constant is given by (atomic units) (26|: 



A = D. e ^ — j . (34) 
Combining equations l32" l l33 l and [34] we obtain that: 

A sc = A|o|. (35) 

The dissipative part of the Liouvillian is identical to equation 1101 (without the tensor product with ly). 

Substitution of equation 1301 and equations 1101 (not in tensor product form) into equation [7] yields a set of equations 
for the matter density matrix elements. In the interaction picture (with Hq given by H m ) and assuming perfect 
matter-field resonance, these equations take the form: 

Poo = 2r i(n i + l)/9n - 2r in iPoo - 2r 02 7io2Poo + 2r 02 (no2 + l)p22 

pu = -iA sc p2i + «A sc pi2 - 2r i(n i + l)pn + 2r in iPoo 

P22 = -«A sc pi2 + «A sc p 2 i - 2T 02 {n 02 + l)p 2 2 + 2r 02 no2Poo 

P12 = -«A sc p 2 2 + «A sc pn - r i(n i + l)pi2 - r 02 (no2 + l)pi2 

pox = «A sc p 2 - r i(2n i + l)poi - r 02 no2Poi 

P02 = «A sc poi - r 02 (2n 2 + l)po2 - r inoiPo2 

fax = P\2 

P10 = PlO 

f>20 = P20- (36) 

B. Thermodynamics of unipartite systems 

Heat flux (Q) and power (P) for unipartite systems with external (time dependent) forcing were originally defined 
by Alicki @: 

Q = Tr(§ff)=Tr{£ d [p]ff} (37) 



at 

(38) 



C. Steady state solution of the semiclassical ED-JCM 



Before we derive the steady state power and heat flux we wish to discuss the main differences between the semiclas- 
sical ED-JCM and the semiclassical theory of the laser due to Lamb [2jJ. Firstly, in Lamb's model the material system 
(the atom) has two levels, whereas in the semiclassical ED-JCM the matter has three levels. Secondly, in Lamb's 
model pumping and decay of the two lasing levels are phenomenological (where the pumping function affects the 
field and thus the interaction term in the Hamiltonian) , whereas in the semiclassical ED-JCM pumping and dumping 
of matter population from the ground state to the two upper lasing levels is achieved through the full dissipative 
Lindblad superoperator. Thirdly, in Lamb's model the field is allowed to decay phenomenologically, whereas in the 
semiclassical ED-JCM discussed in this paper the field does not decay. Finally, in Lamb's model, Maxwell's equations 
for the classical field are solved self- consistently with a quantum perturbative solution of the atomic density matrix, 
whereas in the semiclassical ED-JCM discussed here the field is not accounted for directly. As was mentioned previ- 
ously, in the semiclassical model of Geva and Kosloff the field is not accounted for directly as well. Therefore, cavity 
damping is not incorporated, and negative steady state power in the atom signifies an increase in the field's energy. 

The steady state solutions for poi an d P02 is poi = P02 = 0, since 7P01 = 0, P02 = (3poi, and P, 7 > (after applying 
the steady state condition P02 = Poi = 0). Combining the equations for pi 2 and P21 at steady state (pi2 = P21 — 0) 
yields a central equation: 



IP12 1 cos0(r O i(n O i + 1) + T 02 (n 02 + 1)) = 0, 



(39) 
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where <j> is the phase of the p±2 density matrix element. There are now three possible physical solutions. 
A. |p 12 | = 0. This yields: 

1 

Poo = 



l + 2z 

Pn = P22 = zpoa, (40) 



where z = n n °+ 1 = . Note that this corresponds to a very specific choice of parameters. 

2 



B. cos4> = 0, <f> = ^r- This yields a situation where there is no inversion of the atomic levels: 



Pn - P22 = -^( r oi("oi + 1) + r 02 (n 02 + 1)) < 0. (41) 

Moreover, it leads to a positive atomic steady state power which corresponds to attenuation of the electromagnetic 
field. This is outside the scope of the current paper, and will be explored in more detail elsewhere [28l |. 
C. cos0 = 0, <fi — Tj. This yields a situation where there is an inversion of the atomic levels: 

Pn - P22 = ^(Toi(noi + 1) + r 02 (rio2 + 1)) > 0. (42) 

Asc 

The steady state solutions for the pn, p 22 , and |pi 2 density matrix elements is obtained through the solution of the 
following set of equations: 

r i(2n i + 1) -Toinoi -A sc 
-r 02 n 02 -r 02 (2n 02 + 1) A sc I I p 22 I = I -r 02 n 02 I . (43) 



Asc A S c 

The solution of equation l43l can be written as: 



roi(noi + l)+r 2(no2 + l) 






A/F 
B/F iD/F | , (44) 
-iD/F C/F 



where A, B, C, D, F are given in the Appendix II. 



D. Steady state heat flux and power in the semiclassical ED-JCM 

We are now in position to compare the steady state heat fluxes and power of the fully quantum model with the 
analytical solutions of the semiclassical model. Applying Alicki's definitions (eq. [37] and eq. [38| to the semiclassical 
ED-JCM at steady state yields: 



pss 

Qh 



2r iro 2 A^ e (noi - nogjoJ 
F 

2r ir 02 A^ ; (noi - n 02 )(a;i - u>p) 
F 



a ss 2r 01 r 02 A^ e (n 01 - n 02 )(uj 2 - uj ) 

where F — ^(roi, r 02 , n i, n 02 , A sc ) is a positive constant given in the appendix. We note that at steady state 
Tr{C dH{c) [p ss }V} = 0, and thus Q s £ {c) = Tr{C dH(C )[p ss ]H m }. 

Under the condition A ^> T, the reservoir heat fluxes and power for the fully quantum ED-JCM are found numerically 
to be independent of |a| for the range < \a\ < 10 (which corresponds to an initial coherent state ranging from no 
photons at all to 100 photons in the cavity). There are 0.5% deviations for the higher field strength range (where the 
initial number of photons in the cavity is close to 100) due to a slightly rougher truncation of the Fock space. 

The analytical semiclassical hot reservoir heat flux and power in the range 0.1 < \a\ < 10 are practically independent 
of \a\, and agree almost perfectly with the numerical steady state fluxes in the fully quantum model. However, as \a\ 
decreases below 0.1 the semiclassical reservoir heat fluxes and power change dramatically. This is of course expected, 
as A sc oc \a\, and thus when the field's amplitude decreases below 0.1, X sc is no longer much bigger than T. Under the 
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condition A > T, we find essentially perfect agreement between the numerical steady state fluxes of the fully quantum 
ED-JCM and the analytical steady state fluxes of the semiclassical ED-JCM. Therefore we can state that as far as 
thcrmodynamical fluxes are considered, the semiclassical ED-JCM captures the true physical picture. One important 
exception is that in the semiclassical treatment, if there is no initial field present at all (Eg = 0), amplification can 
not take place. 

For completeness, several thermodynamic 

modes of operation ower: 



(46) 



where = 
is a positive constai 
go together hand ii 
amplification of the 

■ >i 3c r 0i r 02("01-"02) 

E 

Therefore, \p\\\ is a 
plotted in Fig. [6] fo 
more detail in a fort 



□ □ 



□ 



ly when noi > rto2 (E 
■sion and amplification 
■ (corresponding to an 
B reveals that |pf| = 

Mily when no2 > ^oi- 
rence at steady state is 
on will be described in 
) is observed only when 




FIG. 6: Semiclassical atomic coherence, a) A » T. b) A = T. c) A « T. 



A w r. 



E. Engine efficiency 



What about the engine's efficiency? In the previous section we mentioned that the (numerical) efficiency of the 
quantum amplifier always matches the ratio obtained from Scovil and Schulz-DuBois's intuitive definition. Before we 
calculate the engine's efficiency, we wish to obtain Carnot's formulation of the second law in differential form. We 
begin with Spohn's entropy production function: 



SS Q„ Q c 



where ^ is the three-level system entropy change, and — j"'^' is the entropy flux from/to the hot (cold) reservoir. 



ic fVi£i +■ In T"flCi_l Oi-frcA ciTGi- £am Dnfrnmr plnonrrn onrl Q £ 

F H(C) 

At steady state ^ = 0. We now wish to rewrite in terms of Qff and P ss . The quantity E = Ti{pH}, which 
measures the energy flux including the atomic-field interaction energy is given by: 



E = Tr { ^h) + Tr { p^- }=Q H + Qc + P- Us) 



At steady state: 



dt V dt 



E» = QS + Qc + P ss = 2T ^ X ^- n ^~^-^. (49 ) 

F 
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The quantity E ss is zero only at perfect atomic-field resonance. However, the quantity E m = Tr{pH m }, which 
measures the energy flux without the atomic- field interaction energy, and was introduced originally in [9[, is zero at 
steady state, as H m does not depend on time. Expanding E m yields: 

E m = Tr{piJ m } = Qhtti + Qcrn + P m , (50) 

where QH{c)m — 1fc{A£ff(C) [p]H m } and P m = — j-Tr{p[H, m , V(t)]} are the alternative definitions for heat flux 
and power introduced in [9|. At steady state: (1) E m — and hence Qc m = — + Pm)i an d (2) since 

Tr{£dH(c)[p ss }V} = 0, Q S H(c)m = Q S H(GY Therefore we can replace Qq in equation l47l with —(Qh(c) + ^m) where 

v ss 2roiro 2 Af c (7ioi - n 02 )(u)i - uj 2 ) , , 

p m = p , (5!) 

and obtain: 

pss rp _ p 

T,= --&<iZ=- -2, (52) 
Qh Th 

which is Carnot's efficiency formula in differential form. We note that equation [52] is always true regardless of a 
resonance condition. Moreover, we wish to emphasize that rewriting Qq in terms of Q s jf and P" for non-resonant 
cases is possible only through the alternative approach to energy flux in unipartite systems discussed in [9( . 

Substitution of P" (which is identical with P ss at perfect resonance) and Qff into the engine's efficiency formula 
at steady state (^ff- = 0) yields: 

. = -5 = ^ = ^ (53) 



UJ,> 



which is identical with the maser's efficiency defined intuitively by Scovil and Schulz-DuBois. Our model offers a 
statistical description for the reservoirs, and it allows us to derive thermodynamic fluxes, which in turn yield Scovil 
and Schulz-DuBois's efficiency formula. 

We note that Geva and Kosloff @ also considered a semiclassical model for a three-level amplifier. The main 
difference between their model and the semiclassical ED-JCM is that in Geva and Kosloff 's model the time dependence 
of the classical field affects the dissipative superoperator. As a result, the steady state efficiency in the model by Geva 
and Kosloff depends on the power of the field, and hence it is generally not the same as in Scovil and Schulz-DuBois's 
intuitive definition. 



F. Steady state inversion ratio 

In their early work [l[ Scovil and Schulz-DuBois asserted that the ground state population (poo) is bigger than the 
populations in the two excited states (pn and ^22)- This is indeed verified in Appendix II. They also asserted that 
the inversion ratio between the two excited levels is given by: 

t = = e k B T H e k B T c — e k B T H e k B T c . (54) 

P22 

This assertion appears to be well motivated physically, since it would seem that at steady state the population ratios 
between the two excited levels and the pumping level should be related by Boltzmann factors. However, it turns out 
that this is not correct. While it is true that the matter reaches a steady state, as seen in both the quantum and 
semiclassical models, the matter-field system as a whole does not reach a steady state, as was seen by solving the fully 
quantum model in this paper. Moreover, there is no a priori requirement of what steady state populations will be 
attained. We will now demonstrate that the ratio between the two excited levels asserted by Scovil and Schulz-DuBois 
is not correct. Substituting the expressions for the reservoirs' temperatures given in equation 1111 into equation 1541 
yields: 

r = e ln(l/n 02 + l) e -ln(l/n 01 + l) = "-01 (^02 + 1) ^ ^ 

n 02 (n i + l)' 

Substituting Hqi — 10,n 2 = 0.1 yields r = 10. In Appendix II we give analytical expressions for all the density 
matrix elements at steady state, from which a closed formula for r may be obtained: 

p 

r = r(r i,r 2,noi,no2, A sc ) = — , (56) 
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where B,C are positive constants given in Appendix II. Substituting noi = 10,no2 = 0.1, Tqi = To2 = T = 0.001 in 
the analytical expression for r yields (similarly to the quantum model) only a marginal inversion ratio between the 
two excited levels, r = {1.01, 1.00001, 1.00000001} for field strengths E = {0.1, 1, 10}, respectively, where equation 
[55] yields r = 10. 

VIII. CONCLUSION 

We have analyzed a fully quantum model in which a three-level material system is coupled to a single quantized 
cavity mode and two thermal photonic reservoirs in a framework of a heat engine. This gives what is arguably the the 
simplest possible quantum model for light amplification. At the same time, it permits a full thermodynamic analysis. 
Unlike previous work [f| , the field is not considered as an external time dependent force acting on the matter, but it 
is an integral part of the quantum system, allowing us to treat both light and matter on equal footing. We solved 
the ED-JCM master equation numerically, and showed that indeed amplification of the selected cavity mode occurs 
even in this simple model. However, initial field coherence is lost, as seen by the radially symmetric Q function for 
t > (2r) _1 . Moreover, we find that the quantized field mode has an entropy content that changes dramatically at 
short times, and increases very slowly for t > (2r) -1 . The matter-field system as a whole never reaches a steady state: 
at t > (2r) _1 the energy in the field continues to increase linearly in time, which can be analyzed thermodynamically 
in terms of power generation from energy in the hot reservoir. The three- level matter system, obtained by performing 
the partial trace of the full system over the field, does reach a steady state as seen by constant steady state energy 
and entropy. 

Another aspect of the quantum treatment that cannot be dealt with at all within the framework of the semiclassical 
ED-JCM is entanglement. We showed that at short times t < (2F)" 1 the matter-field system is entangled, as seen 
by the negative conditional entropies and the negativity of the partially transposed density matrix. However, at 
longer times t > (2I 1 ) -1 we believe that the matter-field system is classically correlated but not entangled, as the 
conditional entropies (which arc almost equal to the partial entropies) and the partially transposed density matrix 
are both positive. 

Based on our previous work on bipartite systems governed by a time independent master equation [q| we were able 
to derive the fundamental laws of thermodynamics. The first law is obtained both for the full matter- field system and 
for the individual (partially traced) subsystems, using thermodynamical fluxes of heat flux and power. The second 
law of thermodynamics in differential form is guaranteed to exist for the full matter-field system through Spohn's [24| 
entropy production function. We define a new entropy production function a m based on matter thermodynamical 
fluxes. Through <j m we show that at steady state, when the main entropy production is due to heat fluxes from/to 
the heat reservoirs, Carnot's efficiency formula is obtained in differential form. 

A strong motivation for this work comes from an early paper by Scovil and Schulz-DuBois [l[ in which they analyze 
a three-level maser as a heat engine. In their work, they intuitively defined the engine's efficiency as the ratio between 
the maser frequency and the pumping frequency. However, they do not connect this efficiency with explicit expressions 
for work and heat, as expected from a thermodynamical analysis of a heat engine. In our quantized field treatment, 
the efficiency formula of Scovil and Schulz-DuBois was found to be in complete agreement with numerical calculations 
based on thermodynamical power and heat fluxes. 

We have also analyzed a semiclassical version of the ED-JCM. We obtained closed analytical expressions for power 
and heat flux at steady state that are in virtually perfect agreement with those obtained numerically for the fully 
quantum ED-JCM. One may conclude from this that as far as steady state thermodynamical fluxes are concerned, 
the semiclassical model is sufficient. Furthermore, from our analytical results for power and heat flux we were able 
to recover Scovil and Schulz-DuBois 's efficiency formula analytically. One of the assertions in the work of Scovil and 
Schulz-DuBois is that the ratio of populations in the two excited levels is given by a product of Boltzmann factors. 
We showed analytically that this last assertion does not hold in general. 

In future work, we intend to explore further the other thermodynamic scenarios implied by the present model, both 
semiclassically and quantum mechanically. Of particular interest is the reversal of the present mode of operation of 
the engine so that it operates as a refrigerator for light. 
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Appendix I: Derivation of the damped JCM master equation 

The master equation for the resonant Jaynes- Cummings model (JCM) with atomic damping in the interaction 
representation is given by: 

p I af =C h [p I af ]+C d [p I af ], (57) 

where p l a j is the combined atom-field density matrix in the interaction picture, and £h[Paf] an d ^-dip^f] are given 
by: 

= -\[V I ,p I af ] = -i[\(*-ai+<T+a),p I af ] 
Zdlpif] =r(n t/l + l)([ --,p^ -+]+[ --p^, CT +])+rn4[ -+,^ /0 --] + [ CT +p^, --]), (58) 

where cr+(a^) and er~ (a) are atomic (field) creation and annihilation operators (er+er - — <j~ er+ = er z ; er z being 
the Pauli z matrix). A, T, n t h are the atomic-field coupling constant, Weiskopf-Wigner decay constant, and the 
number of thermal photons, respectively. Note that since equation [57| is a master equation of a bipartite system, 
all the operators in equation [551 are implicitly tensor products. For example, cr~^&~Paf is shorthand notation for 
(a + <g) l/)(er_ <8> &f)Paf The damped JCM master equation is usually obtained by adding the Hamiltonian part and 
the dissipative part. We note that van Wonderen gave an analytical solution for the atomic density matrix in the 
damped JCM [29| . and later studied the entropic behavior of the atom [3(|. In this appendix, we derive the full JCM 
master equation by applying the weak-coupling, Markovian and Weiskopf-Wigner approximations, and using a set 
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of unitary transformations. The derivation of the dissipative part follows closely the derivation given by Scully and 
Zubairy \D§. 

We start with the full system (atom-field)-bath Hamiltonian in the Schrodinger picture: 

H = H S + H b + V ab , (59) 

where H s , H b , V sb are given by: 

H s = H a + Hf + V a f = h-^-& z + fuofafa + h\(&~ a 1 * + <r+a) 

k 

Vsb = h^2\ k {&-a{ + *+a k ). (60) 

k 

We denote by s the atom-field system, by b the bath which is composed of an infinite number of oscillators where 
the operators of each oscillator are denoted by subscript k, and by X k the atomic-/cth mode coupling constant. The 
hat notation indicates that all operators are implicitly tensor products with the appropriate identity operators. For 
example <x _ = cr _ ® 1/ <g> 1;,, a = l a ® a <g> 1&, and a k = t a <g> 1/ <E> a k . The above Hamiltonian is written under the 
rotating wave approximation (RWA) meaning that only energy conserving terms are considered. Note that only the 
atom is coupled directly to the bath modes. The evolution of the full system-bath is purely Hamiltonian: 

i " 

Psb = Ch = —r[H,Psb]- (61) 

We now move to the system-bath interaction picture (denoted by superscript I): 

l Psbe 

PL = -^[vLpLI (62) 

where: 

Ho = H a + H f + V af +H b 

Vj b = s£Ai[S-(t)4e*-» t +l+(t)o ke -*" < ] I (63) 

k 

where er _ (+)(t) = e^^ st &~ (+)e~^"^ s *. In the derivation of equation W2\ we made use of the identity [H s ,H b ] = 0. 
A perturbation expansion to second order in V" S 5 yields: 

PL = -\[Vl b {t),pl b {0)] - ± J* dt'[V/ b (t), [V/ b (t'),pi b (*')]]■ (64) 

Consider the weak system-bath coupling limit, that is p sb (t) = p s (t) ® Pb(0) + p c , where p c is any correlation 
between the system and bath which fulfills Tr(,{p c } = (this holds for p sb (t) both in the Schrodinger and interaction 
pictures). In this case the atom- field system evolves according to: 

Pi = Tr b {pf b } = -jTn{[VJ b (t),pl(0) ® p b (0)}} ±Tr„ { jf* dt'[V/ b (t), [V a £(f ), p{ ® p b (0)]]| . (65) 
Note that in equation |6"51 p b (0) = p b (0) and p^ = e^-^ st p s e^ st . The explicit form of equation |6"51 is given by: 

pi = -i^\k(4)i*-(t),pime iUkt 

k 

- f dt'Y^\ k \ k , {ft-(t)&+(t>)p T s (t>) -&+{t')pl(t')a-(t)]e^ t -^' t \ala kl ) 

Jo k,k> L 

+ (t')pi(t') &-(t')pi(t')*-{t) &-(t)pl(t')ir-(t') 

+ pi(t')^-(t')&-(t)]e^ t - i ^' t \alal,)}+R.c., (66) 
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where H.c. refers to all terms on the RHS. The bath density matrix is now assumed to be composed of a product of 
oscillatory modes each being in a thermal state, that is: 



fi 

Pb = n pk\ pk =e fa, + iV fc +D K)(nfc| ' ( 6? ) 

k n k 

where is the average number of thermal photons in the fcth mode. With this assumption equation 1661 reduces to: 

Ps = f dtiJ2^{lv-(t)° + (t')P T s (t')} <T+(t')pl(t')*-(t)}]e^~^n k (68) 
Jo k 

+ [a+(t)a-(t')p T B (t') a-WplWa+We-^-^ink + 1)} + H.c. (69) 
The sum over k is now replaced by an integral: 

y f2Tt pTr roc 2 



where V is the quantization volume and Vk is the fcth mode oscillation frequency. Substituting A| = 2 he v ^ >2 cos2 ^ 
(T> is the transition dipole matrix element, and 9 is the angle between T> and the electric field polarization vector), 
and integrating in the Weiskopf-Wigner approximation (extending the lower limit of the integral over V}~ from to 
—oo, and replacing v k — 2^0^ by u) simplifies equation 1691 

Pi = C d [pi] = -Tn th [a-(t)a+(t)p T s (t)-a+(t)p I s (t)a-(t)] 

- T(n th + \)[a+(t)a-(t)pl(t)-a-{t)pl(t)(T+{t)]+B..c., (70) 

where nth = nk (&0 — w / c ) is the average number of thermal photons, and T — c3 is the decay rate. 
We now move to the system Schrodinger picture: 

p s = e-H^Vfes-W* 

p s = -~[H s ,p a ]+e-i H > t £ d [pl]ei H ° t . (71) 

Using the definitions for d-~^(t) and p 1 (after tracing out the bath) it is easily shown that e~ i H,,t Cd[p I ]eTi H! ' t = 
Cd[ps\- Finally, the master equation for the system in the Schrodinger picture is given by: 

Ps = Ch[Ps] + Cd[Pa] 
£h[Ps] = --[H s ,ps] 

C d [ Ps ] = T(n th + l)([(T-,Ps(T + } + [(T-ps,(T + })+rn th ([(T+,Ps(T-} + [(T+Ps,(T-}), (72) 

where we deliberately omitted the superscript S labeling the Schrodinger picture. 
To summarize, we went through the following path: 

Psb -> pi b -> -> p s - 

The first transition takes us from the system-bath Schrodinger picture to the system-bath interaction picture through 
a unitary transformation. Tracing over the bath under the weak coupling, Markovian, and Weiskopf-Wigner approxi- 
mations leads us to the system dissipative interaction picture. Finally, by applying a unitary transformation we move 
to the system Schrodinger picture. 

To complete the analysis we now move to the standard interaction picture which includes both the Hamiltonian 
and the dissipative parts: 



p\ = e i(H a +H f )t pse j-(H a t+H f )t 

Pi = -~[V,pl] + ei( H °+ H f»£ d [p a ]ei( H »+ H f», (73) 

where it can be shown that: e-^ Ha+H ^ t C d [pa\e^ Ha+H ^ t = C d \p*]. EquationES is now identical with equation 
1581 with subscript s replacing subscript af. 
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Appendix II: Density matrix of the semiclassical ED-JCM amplifier at steady state 



The density matrix for the semiclassical ED-JCM operating as an amplifier (pn — p 22 > 0) is given by: 



( 



Poo 
pn i\p\i 
-i\p\i\ P2i 



) 



( 



A/F 
B/F iD/F 
-iD/F C/F 



) 



(74) 



where A, B,C, D, F are given by: 

A = A sc ro2 + A sc i"oi + A sc r 2n 2 + A sc r in i + r 2r 01 + r 2r 01 n 2 + 2r 2r 01 n i + r 2r 01 n 01 
+ 2r 2rQ 1 no2^oi + r 2rQ 1 n 2Woi + r^roi + 2rQ 2 r ino2 + 2rQ 2 r in 2 + rQ 2 r in i 

+ 2ro 2 r in 2noi + r^Toin^noi 
B = A^ c r 02 no2 + A^ c r in i + T 02 Tl 1 n i + T 02 Tl 1 nl 1 + ^r^^noi + ^m^lino2n 2 m 

+ r^roi^-oi + 2r 02 roino2n.oi + ro2roi n o2 n oi 
c = A^ c r 02 no2 + A^ c r in i + r 02 ro 1 n 2 + 2ro 2 ro 1 rao 2 ra i + r 02 ro 1 n 2no 1 
+ rQ 2 r in 2 + r^Toin^ + rQ 2 r in 2noi + rQ 2 roing 2 noi 

D = As c r 2r i(n i - n 02 ) 

F = \ 3 SC T 02 + A^Toi + 3A3 c r 02 n 2 + 3A^r in i + r 02 r^ + 3r 2r^n i + 2r 02 T 2 01 n 2 01 + 5ro 2 r^no 2 noi 

+ 3r 2ro 1 n 02 no 1 + T^roi + 3rQ 2 r in 2 + 2rQ 2 r inQ 2 + 5rQ 2 r in 2noi + 3ro 2 r iriQ 2 noi. (75) 

A,B,C,E are all positive constants, and since A > B,C => poo > Pn,P22- Thus the population in the zeroth 
(pumping) level is always greater then the population in either level |0) or |1). Since D/E = \pi 2 \, a mathematical 
feasible expression is obtained only if rioi > n 02 . 



